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Abstract 

We present an application of the theory of stochastic processes to model and cat- 
egorize non-equilibrium physical phenomena. The concepts of uniformly continuous 
probability measures and modular evolution lead to a systematic hierarchical structure 
for (physical) correlation functions and non-equilibrium thermodynamical potentials. 

It is proposed that macroscopic evolution equations (such as dynamic correlation 
functions) may be obtained from a non-equilibrium thermodynamical description, by 
using the fact that extended thermodynamical potentials belongs to a certain class of 
statistical systems whose probability distribution functions are defined by a stationary 
measure; although a measure which is, in general, different from the equilibrium Gibbs 
measure. These probability measures obey a certain hierarchy on its stochastic evolution 
towards the most probable (stationary) measure. This in turns defines a convergence se- 
quence. We propose a formalism which considers the mesoscopic stage (typical of non- 
local dissipative processes such as the ones described by extended irreversible thermody- 
namics) as being governed by stochastic dynamics due to the effect of non-equilibrium 
fluctuations. Some applications of the formalism are described. 



Scope 

The paper is outlined as follows: Section 1 is a brief introduction to the problem of applying 
measure theoretical tools to the study of many-particle physical systems, also some recent 
developments in the field are mentioned. We sketch how the probability measures approach 
has been applied to equilibrium systems (states). In section 2 we present an extension of 
such method for the case of non-equilibrium systems (processes) by means of the Choquet- 
Meyer Theorem on continuous measures. In this section two propositions (2.1 and 2.3) are 
made in terms of non-equilibrium stochastic measures as how to deal with systems out of 
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local thermodynamic equilibrium within a modeling approach. The formalism dealing with 
Linear Response with Memory, (physical) Time Correlation Functions. Section 3 includes 
some results of the successful application of these methods on the aforementioned problems. 
Section 4 presents some brief concluding remarks. 

1 States and processes in Statistical Physics 

One of the great unfinished tasks of non-equilibrium physics, whether statistical or phe- 
nomenological is that of finding a general, model-independent way of describing out of 
equilibrium states, similar to the Gibbs-Boltzmann formalism for the statistical-kinetic de- 
scription that leads to a complete thermodynamical definition of systems in equilibrium. The 
Gibbs-Boltzmann formalism also applies to systems close to equilibrium by means of the cel- 
ebrated local equilibrium hypothesis. This hypothesis lies at the heart of both linear response 
theory and linear (classic) irreversible thermodynamics. 

The main problem in finding such a general theory of non-equilibrium states lies in the 
fact that, due to their intrinsic dissipative nature, such systems do not have a unique definition 
of entropy for each state [|50l [T3l |5T| , because of the non- volume preserving property of its 
phase space ||231[TT1|271 . Since phase space volume is not preserved it is, in general impossi- 
ble to define a unique global measure Il59l , derive from this measure a probability distribution 
function and characterize from it a non-equilibrium system. 

Nevertheless some work in that direction has been proved useful in recent years. For ex- 
ample, by extending the works of Sinai ||62l and later of Ruelle and Bowen [[591 141 Gallavotti 
and Cohen [f2Tl 1221 wrote down an expression for a non-trivial measure characteristic of a 
stationary non-equilibrium state obeying this so-called SRB statistics. This measure plus the 
so-called chaotic hypothesis make possible to formulate, at least at a formal stage a proposal 
for a general theory of non-equilibrium phenomena. In this way probability measures give a 
new insight to our mathematical and physical understanding of non-equilibrium. 



1.1 Physical description and measurability 

Let us consider a description of a system in terms of its states. For continuous systems the 
configurations are countable subsets X of W such that X n A is finite for every bounded A C 
R". A state should be defined as a probability measure on the set of all such X. We must 
give to each open set AcR" the probability of finding exactly n particles in A and also the 
probability distribution of their positions. For every bounded open set A and for every integer 
n > 0, let p,\ > be a measure on A n C W w . We shall say that the p\ form a system of 
density distributions (SDD) if they satisfy the following conditions 11501 : 



(i) 
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for every Ac A' 

(n + p)\ 



[l\{x\, X2-, • • • , x n ) — j— j — - f dx n+ \ ■ ■ ■ [ dx n+p fi/J~ p (xi, X2, • • • , x n+p )(2) 

f—i nM Ja'\a Ja'\a 



p=0 n\p\ Ja'\a Ja'\a 
The above conditions imply the normalization condition 



/ dxi . . . dx n p n A [xi,x 2 , ■ ■ ■ , x n ) = 1 (3) 

n J 

The (statistical) correlation functions^ corresponding to the density distribution [i\ are 
defined by: 

, X2, ■ ■ ■ j S^nJ A 4 a P (^i ) (4) 



{n +p)\ 
P =o P- 

Here x\, Xz,...,x n E A. Equation |4] could be formally inverted to obtain: 

L ^ ( 1)P 
pi 



I 00 (— 1)p r 

jJlixx, X2, • • • , x n ) = — 2, ; — / dx n +\ . . . dx n + p p(x\, X2, • • • , x n + p ) (5) 

n\ ~n v\ Jap 



Since the series in the right-hand side of equation [5] may be non-convergent, there could 
be a case when the statistical correlation functions may not determine the state of the system. 
Nevertheless, it is usually asummed that the series converge. 

It is possible to associate with every SDD n\ a state on a B* -algebra A that we will con- 
struct as follows: 

For every bounded open set A cR" let /C A be the space of real continuous functions on 
(W) n with support in A n . Let us call /Ca the space of sequences (f n ) n >o where f n E /C A and 
f n — for large enough n. Let us now define K, = Ua ^-a- 

We denote by T the topological sum J2 n >o (R u ) n of disjoint copies of the spaces (R u ) n . 
If / = (f n )n>o E )C, a function Sf on T is defined so that its restriction to (R u ) n is: 

n n 

Sf(x l7 X 2 , . . . , X n ) = ■ ■ ■ Yl f P ( X h ■ ■ ■ X i P ) ( 6 ) 

p>0 ii i p 

Let us now consider functions on T of the form ip(Sfi, . . . Sf q ) where fi ■ ■ ■ f q E JC and 
<p is a bounded continuous complex function on R 9 . With respect to the usual operations on 



'These statistical correlation functions are different from the physical correlation functions that we will talk 
about later. The former are the moments of the SDD or Ensemble (as the physicist call it), correlating partitions 
of the probability measurable space; the latter are correlations between physical quantities which evolve in time. 
Obviously both kinds of correlation functions are related but their links are non-trivial, and in fact make up for 
a lot of the actual battlefield of the statistical physicist. 
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functions and the complex conjugation operation (*), the <p(Sfi, . . . Sf g ) form a commuta- 
tive algebra A with an involution. The closure A of A with respect to the uniform norm is an 
Abelian B*-algebra 1501 . 

Given an SDD fj%, a state p on A is defined as follows. For each <p(Sfi, . . . Sf q ), let A 
be such that fi,...,f q e JCa- We write 

P(<P{Sfl, ■■■ S fq)) = Ys / dx l--- dx n Pl(x 1 , ...X n )x (f{Sf U . . . S f q ) (7) 
n>0 

fi = fi(%i, ■ ■ -Xn) Vz, this definition (equation IT]) is independent of the choice of the 
volume A and extends by continuity to a state on A. If we call T the set of states defined 
as above (equation [7]) from an SDD; the mapping p\ — > p is then one-to-one onto T. It 
could be seen that the translations of are automorphisms of A. The algebra here obtained 
A is not separable, hence a direct decomposition of invariant states into ergodic states is not 
possible. However, it can be shown that this could be attained if by an indirect treatment [f6TTl . 

We are now on position of studying the time evolution of the just defined states (i.e. a 
physical process) as a particular automorphism of the algebra A. The time evolution of the 
system occupying the finite region A G R u , with respect to a given physical interaction 
$, is defined by a one-parameter group of automorphisms of A\ which associates to the 
observable A 04(0)) at time 0, the following observable at time t (A(t)): 

oo ±n 

e it H<s> (A) A e -«ff«(A) = i n -[H*(A), A} n (8) 

n=0 U - 

here [B, A] = A; [B, A] n = [B, [B, A] 71 " 1 ]. 

It has been proved [1571 that in the so-called thermodynamic limit (A — > oo), for a general 
class of physical interactions $ with associated hamiltonian^ H$ the following limit exists: 

lim e iTH * (A) A e - lrH * {A) = T T A (9) 

A^oo 

Hereafter we will call T T the modular automorphism operator and T T A will be called 
the r-modular evolution of A or just the time evolution of A. 

2 Whose specific nature it will not affect the following consideration. 

3 a Hamiltonian is a function which essentialy gives the energy of the system in terms of the positions of the 
particles and the nature of the interactions between them, in this case it is sufficient to know that is a bounded 
function of the .t, 's. 
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1.2 Equilibrium states, Gibbs measures and KMS condition 

Let us consider a quantum dynamical system described by a von Neumann algebra A (rep- 
resenting the system's physical observable quantities) with a one-parameter automorphism 
group (taking into account the time evolution of the observable quantities). In equilibrium 
statistical mechanics the entropy S of a given state < ■ > (labelled by an observable) is 
defined as: 

S = —Tr plnp (10) 

where p = p/Tr p is the so-called normalized density matrix. We are to consider states lying 
on a constant observable energy surface E = < H > with H the system's hamiltonian and 
the trace taken with respect to some measure H (an ergodic measure). For such systems, a 
global equilibrium state is characterizable as the state of maximal entropy within this set (of 
constant energy E) . As one knows the constant energy and the normalization of the probabil- 
ity matrix p (Trp =1) constraints are introduced in the maximization procedure as Lagrange 
multipliers [43]. Doing so one gets the usual Gibbs states p = e _/3W . 



It turns out that by analyzing the meaning of this entropy (making it coincident with the 
equilibrium thermodynamic entropy) one gets two main results. The parameter j3 is propor- 
tional to the inverse temperature of the system and the ergodic measure H is identical to the 
Gibbs measure p [|68l . This equilibrium condition H = p (usually called Gibbs condition) 
has the restriction of being limited to apply to finite systems pTfl. 



Now let us examine the quantum mechanical description of a system with an arbitrary 
number of degrees of freedom under a constraint r/ from the standpoint of a measure theory. 
The expectation value of an observable quantity A in a state (•),, could be characterized by a 
density matrix p v as follows: 

(A) v = ^^ (11) 

It is possible to introduce in connection with this measure a so-called modular Hamilto- 
nian H* as: 

Pn = e~^ H * (12) 

Here j3* is a number introduced for later convenience (e.g. as a parameter of periodicity 
for the trace of density matrices). The modular evolution of a property A is given by an 
action-like similarity mapping ll68l as: 

T T (A) = e lH * T Ae- in * T (13) 



4 In order to avoid the confusion between thermodynamic entropies and probabilistic or informational en- 
tropies we will use term entropy when referring to the former and entropic-measure when referring to the latter. 

5 Of course it is possible to replace this modular evolution by the action of a Liouvillian super-operator or 
propagator. See for example (TJ, also l52l For a modern introduction to modular operators see 
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The cyclical behavior of the trace (i.e. Trp v {r) = Trp*(r + 1(3*)) gives the following 
condition on the modular evolution between two states N, M: 

(T T (N)M) V = (MT T+ ^(N)) V (14) 

This condition is called the KMS condition [|42]|. This Kubo-Martin-Schwinger (KMS) 
condition (proposed by those authors as an auxiliary boundary condition for Green's func- 
tions) can be adopted as a simple characterizing property of equilibrium states in the alge- 
braic formulation (also called C*-statistical mechanical formulation), which makes sense for 
an infinite system (in contrast to Gibbs Ansatz) and hence enables one to study an infinite 
system directly, by-passing the thermodynamic limit. This is a generalization of the theorem 
concerning Laplace transforms, which says that the Laplace transform of a function of energy 
E, which is zero for E < 0, is analytic as a function of the conjugate variable (here, time t) in 
the upper half-t-plane 0. The KMS generalization of this says that, if the energy of a theory 
is positive then the equilibrium state gives a two-point two-time correlation function that is 
analytic in time-difference in a strip in the upper-half plane of width 1 / K*T, where T is the 
temperature, and is periodic with imaginary period 1/ K*T, K* is a measure-specific con- 
stant. The periodicity follows from the cyclicity of the trace. The justification of this KMS 
condition as the characterization of equilibrium states is that KMS states along with ground 
states and ceiling states (corresponding to ±0 temperature) are precisely those states which 
are stable against local perturbations. 

If the system is in a global equilibrium state (-) v =< ■ >, /?* is just the inverse temper- 
ature j3 = \/K B T and the modular hamiltonian H* is the usual mechanical hamiltonian H. 
In this case modular evolution generates time evolution. Stationarity and dynamical stability 
are minimal requirements for a state to be called an equilibrium state; A passivity theorem 
OTI can be rephrased by saying that it is precisely the KMS states which are distinguished 
by the second law of thermodynamics in Kelvin's formulation: there are no cyclic processes 
converting heat into mechanical work if the state of the systems obeys the KMS condition 

ma. 

Local equilibrium could also be defined through a KMS-type condition. We can think of 
local thermodynamic equilibrium as a state that cannot be distinguished from a global equi- 
librium state by infinitesimally localized measurements 11361 . The quantitative description of 
the term infinitesimally localized measurements could be given in terms of a one-parameter 
scaling procedure ll36l over the observable. We introduce a parametric diffeomorphism as a 
scaling procedure as follows. 



6 For a proof of the theorem in the above context of modular evolution see reference l57l . 
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Definition 1.3 

Uniparametric state space scaling 

A diffeomorphism of the state space that preserves the topological and ergodic character of 
such space and projects a set of points (called quasi-local points) into equilibrium points by 
means of a change of scale is called a Uniparametric State Space Scaling (USSS). 

Let us consider an observable depending on n-points in real space A = A(xi,x 2 , ■ ■ ■ , x n ). 
For the sake of simplicity we will consider a unique field <p spanned by the various values of 
the observable and depending on spacial localization as follows: A = 0(xi)0(x 2 ) . . . 4>{x n ). 
A ^-scaling diffeomorphism (USSS) is given as: 

A 5 A = S(O n 0(x^!)...0(x^ n ) (15) 

with AgA the value of A after a scaling by a factor £,£(£) a scaling function (postulated 
in order for the scaling to be well-defined) and 

X((xi) = x* + £(xi -x*) (16) 

x* is the point to which quasi-local points Xj shrink after localization, x* is called an 
equilibrium point-cell in the sense of Onsager. It is possible to see that xt=o(xi) = x* V? 
and that Xe=i{ x i) — x i so that if the system scales with £ — > we could talk about local 
equilibrium (i.e. all points shrink to an equilibrium one) and if £ — > 1 we are in a highly 
delocalized stage and hence local equilibrium is not attained. Of course repeated application 
of the scaling Ag would define a coarse graining operator in the sense evoked by Ehrenfest 
and formalized by Zwanzig [TTOll . 

It is worth noticing that this scaling procedure generates, for each application, a sequence 
of probability distributions po, pi, . . . p n . . .; here po is the initial random distribution, gener- 
ally taken as a uniform distribution, a fact called Stosszahlansatz (by Boltzmann), molecular 
chaos hypothesis or Initial Random Phase Approximation. Since the USSS procedure is, 
by construction, convergent to the state of thermodynamic equilibrium, should be the 
equilibrium distribution function, i.e. the equilibrium Gibbs distribution, = p. A weak 
convergence condition is generally assumed: 

lim p n = p (17) 

n— >oo 

Even if the state of equilibrium (global or local) is not attained the description of the sys- 
tem under a modular evolution could be done. In this case the measure will depend on the 
space-time localization described by the scaling A^. For a variety of conditions a systemati- 
zation for the scaling effect could be done in terms of a time-dependent distribution function 

7 We call quasi-local to those points with a value of the property A sufficiently close to its equilibrium value 
A* so that after a finite scaling they will indeed take the value A*. Hence quasi-local points are good candidates 
for representing local equilibrium states. 
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fi T (£) for the values of £. The rime r can be considered as the indicator for the direction of 
convergence. In some cases the distribution function is given in the form of a stationary abso- 
lutely continuous uniform measure (or a family of Cauchy convergent continuous measures) 
and so the weak convergence argument -equation 1T7T applies. B4~lfT0"ll . 

Several physically relevant problems could be better understood in terms of uniformly 
continuous measures of parameter state spaces. The problem of extending the notion of 
Gibbs states to quasi-local interactions [45], coupled quantum systems 11371 . the hierarchi- 
cal structure of physical descriptions in terms of the coarse-graining OOll exemplifies such 
problems. The recent arise of the so-called non-extensive statistical mechanics [|64l and the 
somehow dubious q-based entropies f72l are an indirect consequence of a problem that has 
been called non-canonical averaging (3). This problem has been related with the issue of 
aging in physical systems in the form of anomalous phenomena such as the ones surrounding 
the glass transition, a problem that could also be treated as the effect of a non-local in time 
measure. Quasi-local interactions, coarse-graining- scale effects and non-canonicity are all in 
a class of phenomena for which the description in terms of absolutely continuous measures 
(phase space densities) results enlightening. 

As we already noticed a central issue in non-equilibrium physics is the description of ir- 
reversible processes in a formalism akin to the Gibbs -Boltzmann Ansatz. Since we have seen 
that a measure-theoretical description carries on the double duty of characterizing equilib- 
rium (and also local equilibrium) states through a KMS condition, and also providing a tool 
for calculation of the evolution of equilibrium quantities (i.e. modular evolution via an equi- 
librium distribution called the Gibbs measure p) it is natural to think of a possible extension 
of this procedure for the characterization of non-equilibrium processes. Nevertheless, even 
if this goal has been pursued by some of the most brilliant minds in thermal physics: Boltz- 
mann, Landau, Onsager, Zwanzig, Green and many others; the problem has been proved to 
be very difficult to tackle. The reason is that since we have to base our measure theoretical 
description on a non- stationary hamiltonian (a strongly time-dependent quantity), the explicit 
form of the associated measure is given in terms of the solution (impossible in practice) of 
the complete many-body problem. Several approaches have been attempted, ranging from 
projection operator techniques 117011711 to memory functions [3], variational principles [1681 , 
etc. 

In the next section we will introduce an alternative based on a kinetic-theory-founded lfT5l 
extended irreversible thermodynamical formalism 11131 171 fT6l taking into account the fact that 
in the typical time and length scales of extended non-equilibrium phenomena (the so-called 
mesoscopic stage) the fluctuations of the macro-variables will play a very important role. 
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2 Non-equilibrium measures 

By application of the so called modified moment method for the kinetic theory extraction of 
macroscopic non-equilibrium information, it has been proved |[T2ll that a probabilistic mea- 
sure (a "quasi-gibbsian" measure) exists for an arbitrary pair of (non-equilibrium) steady 
states given in terms of the entropy production (more precisely the uncompensated heat pro- 
duction) at the steady states. This measure gives rise to an extended Gibbs relation that 
generalizes the Gibbs relation for the entropy change usually employed in equilibrium ther- 
modynamics 0. 

The distribution function f ne — f + f™ e + fg e + . . ., if examined under the so-called 
functional hypothesis |[T5l [TTll gives rise to an irreversible thermodynamical formalism that, 
in spite of being of a nonlocal and nonlinear (far from equilibrium) nature obeys a form of 
canonical thermal averaging; this fact will be very useful later in this study. Here fo is the 
equilibrium distribution function and f™ e , i = 1, 2, ... are non-equilibrium contributions to 
the distribution function also called higher order moments of the distribution function. 

The specific form of the measure p depends on the solution of the kinetic theoretical 
description of the system under study. Once again, except for a few ideal systems like the 
classical and quantum diluted gas there is no closed solution for the kinetic equations. How- 
ever, since the measure depends on the kinetic solution in terms of a distribution function 
l 1 = l 1 (/o) f\ e i f?, e i ■ • •) ( even without having the explicit solution for it), it is possible to 
see that the passage from the so-called chaotic stage (given by Boltzmann's stosszahlansatz) 
to the fully-developed stationary distribution shall induce a dynamic convergence sequence 
in the measure. Given these facts the main contribution of this work could be summarized in 
two principles: 



Proposition 2.1 

The measure p T associated with the state space modular evolution T T (X) for a system out of 
(local or global) thermodynamic equilibrium and described by variables X(t) is given by a 
sequence of time dependent distribution functions. This means that we are able to recognize 
the measure as the stationary solution of a stochastic process Q T (t). 

More explicitly we can state that: 

lime T (i) = /i r (18) 

T T A= [ e T (t)F(A(t),r) dt= f p T F T (A(t)) dt (19) 

F is called the kernel or modulus of the given modular evolution T and is a still-undefined 
continuous monotonic function of the state space fields X. The dynamic variable A is an 

8 System-dependent in the Physicists saying 
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arbitrary function of the state space fields X, A(t) = A (Jc(tfj, T is a time-continuous 
stochastic process. The sequence induced by equation [191 converges according to the distri- 
bution Q T already mentioned in connection with the scaling procedure in the previous section. 
Although the explicit form for the distribution function for the scaling parameter Q T is a pri- 
ori unknown we will see that is homeomorphic to the convergence sequence for the kinetic 
distribution function f ne . 

Theorem 2.2 

The stationary state of the Stochastic process T corresponds to the non-equilibrium prob- 
ability measure \i at least in the sense of means {i.e. almost-surely) according to equation 

dH. 
Proof 

Let VL T (£) be a distribution of r-time dependent scaling factors in the sense ofUSSS. Let us 
assume that the distribution is stationary and has a compact support (i.e. it is dense in state 
space). By construction of the USSS, the stationary distribution corresponds to an equilib- 
rium (local or global) distribution of state space points. Now let fj, = /x (/o, fi e , Z^ 6 , . . .) be 
an e-dependent or kinetic distribution, where e is called the uniformity parameter or Knud- 
sen 's parameter. By definition p is stationary and has compact support in phase space. Since 
the stationary distribution \x (/o) is an equilibrium (local or global) distribution there exists 
a transformation from one equilibrium distribution to the other (i.e. in the thermodynamic 
limit both distributions are equivalent). 

This condition suffices, to take into account the fact that as £ goes to zero the kinetic 
distribution f ne should converge to the equilibrium distribution f . This is so, since the con- 
dition of equilibrium ( in the sense of spatial homogeneity of the thermodynamic potentials ) 
gets attained (see equationU6\). notwithstanding the difference in the speed of convergence of 
both distributions. 

In some sense, Theorem 2.2 is a kind of informal version of the Monotone-Convergence 
Theorem ll67l . The second proposition, namely equation ( fT9l follows Theorem 2.2 plus the 
assumption of modular evolution. Let us see this in terms of measures: 

Let F be a locally convex topological vector space and K a convex, compact subset of 
F. The dual M. of the closure of K consists on the measures on K. Denoting by A4 + the 
convex cone of positive measures and by M. x the set of positive measures of norm 1 (i.e., the 
probability measures). If p £ Mi there exists p £ K, such that: 



f(p) = J f(a)dp(a) 



(20) 
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p is usually called the resultant of p, §5$. If p, E M.\ has a resultant p then p can be 
approximated weakly by measures p! G M.\ with resultant p and finite support ^ 

If we take this result in terms of equation [T8l equation [19] follows directly. In the second 
expression at the right hand side of equation [19] has been used the so-called tower property of 
conditional probabilities ll67l . 

We have still retained the notation of discrete evolution in order to show up the analogy 
with usual (i.e. finite volume) Gibbs measures, however we must note that p T F T represents 
the operation of a time continuous stochastic process r (i) over a field through the action of 
a r-continuous semigroup F 0. 

Once we have defined this measure we are able to make explicit assumptions to what the 
time evolution of a macro-variable will be in terms of correlation functions under the context 
of non-local irreversible thermodynamics. 

Proposition 2.3 

The time evolution of a (non-local) irreversible thermodynamic field Z given in terms of a 
modular evolution operator Y T generates a form for the 2-time correlation function given by 
the r-continuous stochastic map: 



Similar expressions could be write down for higher order and crossed correlations. 

As it could have been already noticed Propositions 2.1 and 2.3, are an extension of linear 
response theory for the (non-linear) case of a time dependent measure. At the moment and 
having recognized the impossibility to derive this measure from microscopic models we have 
assumed it as an stochastic variable. Of course if this stochastic measure converges (at least 
in the sense of means) -see equation \T7\- to the Gibbs measure p we recover the usual linear 
response theory. 

2.1 Linear response with memory 

Non-equilibrium thermodynamics often addresses the problem of transport through material 
media. Hyperbolic transport equations (such as the MCV (6] [651 equations and the telegra- 
pher's equation j26J) taking into account the lag on the response due to the finite velocity of 

9 A formal proof of this Choquet-Meyer Theorem (too lengthy to be included here) is given in reference 
1() Since F T represents the time evolution of a mesoscopic and hence possibly dissipative system the inverse 
operation may not satisfy existence and/or unicity so we will refer to the effect of operating F as the action of 




(21) 



a semigroup. 
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perturbations have been derived from several points of view. Ranging from phenomenologi- 
cal arguments to purely microscopic transport E6ll and also probabilistic methods such as the 
persistent random walk |[2~8ll . On the other hand linear response theory represented a pow- 
erful theoretical tool to cope with transport phenomena from a dynamic and thermodynamic 
standpoint. In the last years linear response functions incorporating memory (i.e. the effect 
of the lag in the response to an applied field) appeared in such problems as delayed transport 
in electronic devices, molecular hydrodynamics and rheology of structured fluids, and also 
solid state phenomena such as the dynamics of Abrikosov vortices. 

— * 

Let us first consider two related dynamic variables say the magnetic field ri{t), and the 
magnetization M.. At constant temperature in the low-field limit there is a linear relation 
between the magnetic field and the magnetization: 



Here Xt 1S a tensor response coefficient called the susceptibility tensor. However, it 
is known that it takes some time t m for a material media to achieve magnetization under 
exposure to a magnetic field. In other words, there is a lag in the response to the field. This 
effect is due to the fact that magnetization is a non-equilibrium process. If we look at the 
effect of mesoscopic fluctuations of the magnetic field due to non-equilibrium evolution in 
the sense already evoked we will see a possible explanation. Since the magnetic field is 
in disequilibrium we can assume it undergoes a stochastic time evolution. It is possible to 
associate an automorphism T v on the magnetic field with this stochastic evolution: 



But since the modular automorphism is induced by a time-continuous measure fi T we 
have: 



By proposing a form for the stochastic process associated with the measure fi T of the 
automorphism T v we can model this non-equilibrium process. If we look for the behavior in 
a weak stochastic limit we are able to consider a first order stochastic process. Let us consider 
r,, as a modular automorphism whose associated stochastic process is a simple exponential 
decay in time since this is a prototypic first order stochastic process ||39l . In this case equation 
[25] reads: 



M = x°t H(t) 



(22) 





(25) 




(26) 
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here X P is the amplitude of the associated distribution and t p is a characteristic time 
(akin to Poisson's time). If we define a non-equilibrium susceptibility Xr e ( U P to first order) 
as Xr e = Xt^p> an ^ recognize Poisson's time as the relaxation time associated with the 
non-equilibrium process of magnetization we get: 



Equation [27] is the usual expression of a linear response with exponential memory (also 
called fading memory). If we look at the differential representation of equation [27] we will 
find the hyperbolic MCV-type transport equation. It is worth noticing that even if we start 
with a linear response relation and the measure is induced by a first order stochastic process, 
the resulting equation possess a non-markovian character, i.e. memory in the response. If 
the time scale of the system is much more slower than its stochastic time (t c h ar >> Tp) 
then we can consider the limit Tp — > 0. In this case the exponential decay distribution is 
just a (^-distribution and the usual linear response without memory is recovered. By means 
of the analysis of the associated measure (and its related stochastic process) one is enable to 
perform a selection of the effects that we will take into account in a very transparent way. This 
will be a very important issue when considering a collection of irreversible processes taking 
place in a non-equilibrium system, since in this case the presence of irreversible couplings (a 
question related to the relative relaxation times of the various processes) could be cope-with 
in a systematic way. 

2.2 Time correlation functions 

Time correlation functions are usually considered to represent the mean behavior of a large 
set of microscopic dynamic variables under certain averaging assumptions (canonicity, num- 
ber density conservation, energy density conservation, etc.) usually taken into account by 
some kind of Lagrange multiplier formalism (in the case of extended thermodynamics these 
multipliers are introduced, for example, by Liu's Method) [fT9ll38ll44ll . Nevertheless by ap- 
plying a novel irreversible thermodynamics formalism [14] it has been stated that they also 
represent Gibbsian ensemble averages of a collection of macroscopic field variables, when 
the latter are considered (as is the case in extended irreversible thermodynamics [fP3~ll ) as 
coming from stochastic processes (in general non-markovian noises) in a mesoscopic length- 
scale lfT8l[T2l . This statement has only been proved formally for dilute gaseous systems by 
means of quantum kinetic theory |fT3llT8l[T2ll . but there is strong evidence supporting its va- 
lidity for several other physically significant systems. We will use this statement as an ansatz. 

In order for the field averages to be of a more general character, we will use the equiv- 
alence between 2-time correlation functions and propagators 117011711 . or time evolution op- 
erators acting on a dynamical variable. The dynamic variables are the set of extended ther- 
modynamical potentials over its gibbsian set. Time evolution operators are obtained as inner 
products with dual vectors of the aforementioned thermodynamical potentials. From the the- 
ory of stochastic processes Il39ll48l we obtain the statistical properties of such inner products 




(27) 
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by introducing a probability measure ji (or weighting function) in the last stages of study of 
the problem. 



The 2-time auto-correlation function for a dynamical variable A will then be given as (see 
eq.EB: 



As we already stated we are in a position to see the measure ji T as a dual complement of 
A, i.e. (i T A = A^ with (A(t), A(t)*) = \\A(t)\\ 2 J63l. We could also look in a very similar 
way to 2-time crossed-correlation functions and also to higher order correlations. Let us look 
at an example coming from extended irreversible thermodynamics 041 . 

According with the formalism of Extended Irreversible Thermodynamics (EIT) [ 38 ] the 
time evolution of the (entropy-like) compensation function $ is given by: 



We see that equation (|29l is nothing but the formal extension of the celebrated Gibbs equation 
of equilibrium thermodynamics for the case of a multi-component non-equilibrium system. 
This extended Gibbs relation is brought about in the theory by imposition of some consis- 
tency conditions on the non-equilibrium part of the distribution function 031 |T8l |T5]|. The 
quantities appearing therein are the standard ones, T is the local temperature, P and V the 
pressure and volume, is the chemical potential for the species "z", etc. X, and $j are 
extended thermodynamical fluxes and forces. 

We will look up for the effect that a non-equilibrium process, say mass flux will have on 
the thermodynamic description of the system. In the case of a binary fluid mixture we will 
take our set of relevant variables Q — S U T to consist in the temperature T(r, t) and con- 
centration of one of the species C^r, t) fields as the slow varying (classical) parameters set S 
and the mass flux of the same species J 2 (r, t) as a fast variable on the extended set T . For the 
fast dynamic variables (such as the mass flux J 2 and its conjugated thermodynamic force X) 
characteristic times are much smaller than for the so-called conserved fields (mass, energy 
and momentum densities, etc.) so, the effect of fluctuations in this variables will be greater. 
We will take this fact into account by associating to this fields a memory in the response akin 
to the memory described in the last section. 

We propose linear constitutive equations with exponential memory kernel (see section 
2.1) for the following reasons: a) The associated transport equations are hyperbolic (of the 
Maxwell-Cattaneo- Vernotte type) (6l |65l so causality is taken into account, b) These hyper- 
bolic transport equations are compatible with the postulates of EIT Il38l 1441 [531 [541 1551I . c) 
Similar equations can be derived for coupled non-markovian stochastic processes |[66ll47l and 




(28) 



Td t m = d t U + Pd t V - £ Tidtd + dt$ 



(29) 
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since stochasticity has been associated with the major role of fluctuations in the mesoscopic 
description level of EIT the outcome will improve our understanding of phenomena occurring 
in such mesoscopic stages, d) The stochastic process associated with this set of two coupled 
constitutive equations called a semimartingale has been extensively studied and its known to 
accept absolutely continuous measures a fact that will become very useful lTT0l 1671 1331 . 

The constitutive equations are therefore chosen to be, 



The Aj's are time-independent, but possibly anisotropic amplitudes, u is a unit vector 
in the direction of mass flow and r/s are the associated relaxation times considered path- 
independent scalars. Since we have a linear (thought non-markovian) relation between ther- 
modynamic fluxes and forces some features of the Onsager-Casimir formalism will still hold. 
The resulting thermodynamic potentials and its derivatives are assumed to belong to the 
quasi-gibbsian set of time dependent macroscopic functions mentioned in reference lfT2l . 
over which we will perform a thermodynamic average (i.e. a phase field average as mea- 
sured by a quasi-gibbsian measure characteristic of the non-equilibrium thermodynamic state 
space) in order to obtain the 2-time dependent correlation function for, say the concentration- 
concentration correlation function (Ci(f, t), Ci(r, £')). If the thermodynamic average is per- 
formed under isotropic canonical conditions Il68l 1691 1321 the resulting correlation function is 
given by: 




(30) 




(3D 



(C 1 (r,t),C 1 (r,t>))= I d(f,t) ^jdff.iVi' 



(32) 



For the case being treated it was showed [34] that equation[32]is given as follows: 



(C 1 (f,t),C 1 (f,t'))= [ S(r) e-"(*+*') W{t-t')dt' 



(33) 



with 



E(r) = AB + B 2 



(34) 




(35) 




A 1 A 2 r 1 iJ(r)T(r) 
n(l — mi) 2 



][Air 1 2 -r 2 ( 



(1 - nn) 



AiTi 



+ 



Ai(l-nn)-l 
(1 - nr 2 ) 



)] 



(36) 



W(t — t') is a stochastic process representative of the measure fift-tf)'> hence the presence 
of non-equilibrium correlations is taken into account by means of this stochastic evolution of 
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the measure0. Equation (1331 is the non-regular part of the composition field time correlation 
function (i.e. the value of the correlation function near the critical point after scaling of the 
thermal and concentration fields). C p is the heat capacity, T(r) and H{f) are known contin- 
uous functions for the amplitudes of the temperature field and the chemical potential field. 

For a broad family of stochastic measures equation (1331) asymptotically converges to a 
limit given by ll29ll : 

(d^O.Ci^O) =«S(f) e~ 2nt (37) 

where k is a constant depending on the explicit stochastic measure under considera- 
tion. For an unitary Ornstein-Uhlenbeck measure (W{t - f) = e 1 *-*' 1 ) convergence im- 
plies k = As we shall see this distribution resulted very appropriate to model this 
kind of non-equilibrium critical system [|34l . In the case of an unitary Gaussian measure 
(W(t — t') = e^~* ) ) for example, an asymptotic solution could not be given in terms of an 
exponential decay, a typical feature of the critical decay of fluctuations. In this Gaussian case 
the solution consist of two contributions: an exponential decay plus an error function type 
mode. This difference eliminates the short-time plateau (known as critical slowing down of 
fluctuations) present in correlation functions obtained under colored noise measures. 

Equationl37lis a form of the van Hove expression for the density-density correlation func- 
tion, also called a dynamic structure factor. 

One of the advantages of this formalism with respect to others, such as mode-coupling 
theory is that we can test different kinds of stochastic relaxational couplings (delta-correlated, 
quasi-markovian, gaussian, lorentzian, Ornstein-Uhlenbeck type, non-markovian, simultane- 
ous non-linear multi-modal coupling, etc.) just by changing the weighting functions W(t—t') 
(that is the family of probability measures jj, T ), within a unifying thermodynamic scheme. 



2.3 Hysteresis in non-equilibrium systems 

It is possible to describe how the phenomena of hysteresis could be well understood in the 
context of an extended irreversible thermodynamic formalism by pointing out that hysteretic 
phenomena is a consequence of the existence of dissipative internal processes in the system 
as was showed elsewhere 11351 . The different contributions to the uncompensated heat pro- 
duction can be explicitly incorporated in the description by characterizing each one by its 
relaxation time. In order to do so it is important to take into account the effect of dynamical 
coupling. One way of taking this effects into account is by examining the coupling of the 
associated stochastic measures, then looking for processes whose measures have stochastic 

"in the case studied then (reference ll34l ) the convergence sequence was Cauchy, with stationary Ornstein- 
Uhlenbeck character for linear non-markovian statistics except in the case where stochastic resonance was 
present. 
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times of the same order, since in this case a dynamic coupling is possible. By enslaving those 
stochastic processes to the faster one it is possible to develop a hierarchy of measures. In 
this case the analysis of relaxation times is also physically clearer than by using the standard 
hysteresis operators |[56l @9]| 

3 Some successful examples of the results of the application 
of the formalism 

Since we introduced the non-equilibrium stochastic measure ji T in Propositions 1 and 2 as 
a tool for modeling systems out of local or global equilibrium its physical validity will be 
ultimately conditioned by the results given by its use against experimental data. We will give 
here some brief results of some cases where its application has led to satisfactory results. Let 
us first consider the extended thermodynamical description of criticality as given in ref. 11341 . 
If we compare the results for the temperature dependence of the diffusion as obtained from the 
calculated 2-time concentration autocorrelation function in the neighborhood of the critical 
consolute point, as given by equation [331 under a suitable measure with the experimental light 
scattering spectra as reported and by looking at the results there, it is clear that the stochastic 
measure was a reliable tool for the description of a highly complex non-equilibrium property. 

For the case of non-equilibrium couplings and its relation to the observed phenomena of 
hysteresis (cf. section 3.3) it has been possible to show ll35l that if the relaxation times of 
the phenomena are compared to the stochastic time given by the associated measure, then a 
criteria is given as to when will we experimentally observe hysteretic phenomena. 

In brief; by applying simple and reasonable modeling approximations within the formal- 
ism outlined in the preceeding sections, we were able to reconstruct the hyperbolic transport 
MCV equations (cf eq. [2Tb of linear response theory with memory. Also a well known re- 
sult of condensed matter theory and phase transitions was obtained, namely the Van Hove 
equation for the density autocorrelation function (eqJjTT). Finally some previous results on 
the phenomena of hysteresis in terms of relaxation times resulted firmly grounded within the 
framework presented in this paper. 

4 Concluding remarks 

In equilibrium statistical mechanics a measure theoretical representation (Gibbs formalism) 
has provided a great insight and also powerful tools for the characterization of equilibrium 
states (via Gibbs or KMS conditions) and the calculation of equilibrium quantities (via ther- 
modynamic averages, i.e. averages given by equilibrium measures). Similar tools has been 
derived for systems close to equilibrium (in local equilibrium) by means of linear irreversible 
thermodynamics and linear response theory. Although this task has been impossible to take 
further away from equilibrium, by using a similar approach we have been able to model 
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expressions representing some non-equilibrium situations in a non-local regime. The appli- 
cation of this formalism (based on the assumption of stochastic evolution at the mesoscopic 
level of description) permits a systematic study of such dissipative systems and in the cases 
outlined give rise to good agreement with available experimental data. Much work has to be 
done in order to fully understand the consequences of this and related formalisms, such as 
the ergodic dynamic systems approach of Gallavotti and Cohen []2T| and the (much debated) 
fractional calculus approach of Tsallis ll64ll . Acknowledgements: One of the authors (E.H.L.) 
greatly appreciates the fruitful discussion with Ruben D. Zarate on modular automorphisms 
in arbitrary spaces. 
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